Spin filtering due to quantum interference in periodic mesoscopic networks 
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We present several new results, extending our recent proposal of a spin filter based on a tight- 
binding model for a periodic chain of diamond-like loops [Phys. Rev. B 78, 125328 (2008)]. In 
this filter, the Rashba spin-orbit interaction (which can be tuned by a perpendicular gate voltage) 
and the Aharonov-Bohm flux (due to a perpendicular magnetic field) combine to select only one 
propagating ballistic mode. For this mode, the electronic spins are fully polarized along a direction 
that can be controlled by the electric and magnetic fields and by the electron energy. All the other 
modes are evanescent. Generalizing the square diamonds into rhombi with arbitrary opening angles, 
we find that increasing these angles widens the parameter range for efficient filtering. A different 
gate voltage on the two sides of each rhombus is found to delocalize the electrons for energies on one 
side of the band center. We also compare our tight-binding model with models which use continuous 
quantum networks of one-dimensional wires, and find coincidence only when one chooses particular 
site energies at the nodes of the network. 



PACS numbers: 

I. INTRODUCTION 

Future device technology and quantum information 
processing may be based on spintronics^, where one ma- 
nipulates the electron's spin (and not only its charge). 
Here we address attempts to build mesoscopic spin fil- 
ters (or spin valves), which generate a tunable spin- 
polarized current out of unpolarized electron sources. 
Much recent effort in this direction uses narrow-gap semi- 
conductor heterostructures, where the spins are subject 
to the Rashba3 spin-orbit interaction (SOI): in a two- 
dimensional electron gas confined by an asymmetric po- 
tential well, the strength of this SOI can be varied by 
an electric field perpendicular to the plane in which 
the electrons move^. An early proposal of a spin field- 
effect transistor^ used the Rashba SOI to control the spin 
precession of electrons moving in quasi-one-dimensional 
wires. 

Some of the most striking quantum effects arise due 
to interference, which is best demonstrated in quan- 
tum networks containing loops. Indeed, interference 
due to the Rashba SOI has been measured on a 
nanolithographically-defined square loop array^. Re- 
cently, several theoretical groups proposed spin filters 
based on a single loop, subject to both an electric and a 
magnetic (Aharonov-Bohm (AB)^) perpendicular fields 
(e.g Refsiii2ii2iiiS). However, such devices produce a full 
polarization of the outgoing electrons only for special 
values of the two fields. Later work considered the ef- 
fects of the Rashba SOI on the conductance of chains of 
loops. These included studies of chains of diamond-like 
loops^^'^^, and studies of finite chains of circular loops^-. 
Although both studies showed some destructive interfer- 
ence due to the SOI, they did not concentrate on the 
tuning of the fully polarized spins which can be obtained 



in certain parameter ranges. 

RecentljJ^, we proposed a spin filter based on a pe- 
riodic chain of diamond-like square loops, connected to 
each other at opposite corners [see Fig. [T]. Unlike the 
above earlier papers, which used a continuum descrip- 
tion for the wires on the network, we used a simple tight- 
binding model, with quantum dots (or 'atoms') only at 
the nodes of the square diamonds. This allowed us to 
obtain transparent analytical expressions for the ballistic 
conductance through the chain and for the outgoing spin 
polarization. We found that a combination of both the 
Rashba SOI and the AB flux through each loop can result 
in destructive interference, which can block the transmis- 
sion of all the spin components except one, which is then 
polarized at a tunable direction. Technically, this sin- 
gle spin direction is associated with a single propagating 
wave solution of the Schrodinger equation, while all the 
other solutions involve evanescent modes. 

Here we extend our analysis of this diamond-like chain 
in several directions. First, we replace the previous 
square loops by rhombi, with a general angle 2(3 (Fig. [T]). 
It turns out that the filter is more efficient for (3 > tt/A. 
Second, we generalize our previous study, by allowing dif- 
ferent site energies (controlled by appropriate gate volt- 
ages) on the various sites in the unit cell (i.e. sites a, b 
and c in Fig. [T|). Different site energies at sites b and c 
turn out to have drastic effects on the ballistic conduc- 
tance. Third, we propose using this filter at fixed electric 
and magnetic fields, controlling the outgoing polarization 
using a gate voltage. (In Ref.^^ we worked at fixed en- 
ergy, and varied the magnetic and electric fields.) Fourth, 
we replace each edge of each rhombus by a tight-binding 
chain of 'atoms' (or quantum dots). In this context, we 
compare our tight-binding approach with the continuous 
quantum network approach used in earlier work on the 
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same geometr^ r^^i^^ . As we discuss elsewherei^, the two 
approaches are not equivalent. 
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FIG. 1: Chain of diamonds. 

Section 2 outhnes the tight-binding model which we 
use for solving the Schrodinger equation on the periodic 
chain of generalized diamonds. Section 3 presents results 
for the polarization of the electrons in the regions where 
they are fully polarized. We compare our tight-binding 
approach to earlier continuum network models in Sec. 4, 
and summarize our results in Sec. 5. 



II. GENERALIZED TIGHT-BINDING MODEL 

The basic theory was presented in Refi^. However, 
for the generalizations introduced here we find it eas- 
ier to use the (rotated) coordinate axes shown in Fig. 
[TJ Setting the opening angle of each diamond to 2/3, 
the coordinates of the sites in the n'th unit cell become 
r(a„) = {nL, 0, 0), r(&„) = r(a„) -I- (L cos/3, L sin/3, 0) 
and r(c„) = r(a„) -I- (Lcos/3, — Lsin/3,0), where L = 
2Lcos(/3) is the basic step along the x— axis. We start 
with the simplest tight-binding model, which has 'atomic' 
sites only at the corners of the diamonds (this will be ex- 
tended below) . The hopping 2x2 unitary matrices within 
the n'th diamond thus become 

[/^^,(n) = e*("+i)'^/2+'"'^i, (1) 

where a' denotes the site Qn+i, cr is the vector of Pauli 
matrices, a = kgpL {kgQ represents the strength of the 
Rashba term) , — sin (3 — Oy cos /3, (T2 = ^i^i P + 
Uy cos /3 and (/'/(27r) = BL^ sin(2/3)/$o is the number of 
flux units through each diamond. For each bond, the 
hopping matrix must be multiplied by the hopping en- 
ergy, Juv — J (below we measure all energies in units of 
J)- 

Within the n'th diamond, the Schrodinger equations 
for the spinors tpain), tphin) and tpc{n) are 

(e - eJV^Jn) - -j{UMMn) + U,,(n)Mn) 

+ UlAn ' l)Mn - 1) + Ul,{n ~ l)^,(n - 1)), 

(e - e,)V'b(n) = - j(C/],(n)^,(n) + (n)V'a(n + 1)), 

(e - eJVc(n) = -J{Ul{n)Mn) + U,A^)Mn + !))• 

(2) 



Except for the special energies e = e^, e^,, which represent 
dispersionless solutions (not shown in the figures), we 
express ipbin) and ipdn) in terms of ipa{n) and ipa{n+l), 
and substitute into the equations for ipain). We end up 
with effective one-dimensional equations, 

4A*a(n)=Wt^,(n-l)+W*,(n+l), (3) 

with 4A = e - e„ - ~ 2j^, 7^. = JV(e - e^), j = b, c, 
and with the non-unitary 2x2 matrix 

EE 2{d-iby(Ty - fe^crj; 2d ^ a^[c^ - s2cos(2/3)], 
2by = 2a^cs cos /3, 2b^ = -~ia_s^ sin(2/3), (4) 

with c = cos a, s — sin a and = 7^6^"^/^ ± 7^6**^/^. 

Assuming a propagating wave, '^a{n) — Ce*^^"x((7), 
where xil) is a normalized spinor, we find that this spinor 
must obey Hx — Ax, with 

77 = (e-^«-^W^ + e^9^W)/4EE (A + B-cr), (5) 

where A — cos{qL)did ~ sm{qL)^d, — 0, By — 
cos{qL)'^by + sin{qL)diby, and B^ = sm{qL)'^bz — 
cos{qL)dibz- Therefore, x is an eigenstate of n • cr, with 
the unit vector n — B/|B|: n • cr^^ = MX^, M = 
and we have A = A + /.t|B|. Thus, e and q must obey 
the equation (A — A)^ = By + B^. At fixed q, this is a 
polynomial of degree 6 in e, so that one expects 6 energy 
bands. In the special symmetric case where e^, — P^'^" 
tially treated in Ref J^, two of these solutions represent 
dispersionless solutions at e = e;,, and thus one is left 
with only four bands (whose shape and location varies 
with efc). The value of only represents a shift in en- 
ergy, so we fix = 0. For > = the dispersionless 
modes become dispersive, and in general the spectrum 
also becomes asymmetric with respect to g <-> —q and to 
e — e. An example of this phenomenon is shown in 
Fig. [21 for the square diamond (/3 — n/A). 

As explained in^**, the ballistic conductance at a given 
energy e is equal to G = {e^ /h)g, where g is the number of 
propagating wave solutions which move in one direction. 
To study (7, we now fix e and solve the spectrum equation 
for q. In the general case, this equation turns into a 
quartic equation in cos(gL). This equation reduces to a 
quadratic equation for = e^,. Out of the solutions for q 
we count only the right-moving propagating modes, with 
real q and with a positive velocity v ~ dt/ dq. As iiiii, 
we again find ranges of e with g — 0, 1 or 2. Generally 
g depends on all three parameters e, a and 4>. Unlike 
Ref.^**, where we fixed e and presented results as functions 
of a or (j), here we fix 4> and Fig. [3] presents contour plots 
of g{e,a/'K) at (/) = 7r/2. As already indicated by Fig. 
[21 changing can open the large gap which existed in 
the symmetric case near e = 0. This can be seen by 
comparing the top two plots in Fig. [3l the right hand 
side plot, for e^^ = J and — 0, is asymmetric with 
respect to changing the sign of e, and it exhibits non- 
zero ballistic conductance at small positive energy. The 



gate voltage governing can thus be used efficiently to 
vary the ballistic conductance between zero and non-zero 
values. 




FIG. 2: Spectra for /3 = 7r/4, </> = An, a = .2n. Four full 
lines: e,, = = 0. Six dashed lines: ej, = .5, = 0. All 
energies are in units of J. 









III. SPIN FILTER 

In the regions with 5 = 1 we have only one propagat- 
ing mode. As explained above, each mode is associated 
with a single spinor xil)- The spinor of an electron com- 
ing from the left will in principle be written as a linear 
combination of all the four (or six) solutions. However, 
when g = I all the right moving modes except one are 
evanescent, i.e. they decay with distance. Therefore, 
the regions with g ~ 1 represent full polarization of the 
conducting electrons. For the symmetric case. Fig. [3] 
also compares between three values of the rhombus angle 
f3. Interestingly, the square diamond (/3 = 7r/4) is not 
the optimal filter; increasing f3 broadens the regions with 
g — 1, where we have full polarization of the electrons. 
Therefore, we present below results for /3 = .35tt. 



FIG. 3: (Supplied separately) Contour plots in the e — a/n 
plane of the ballistic conductance g. (a) /3 — 7r/4 and £(, = 0, 
(b) /3 = tt/4 and = J. (c) ej, = and /? = .IStt. (d) ej, = 
and P = .357r. All other site energies are zero, and all plots 
have (j) = 7r/2. The values g = 0, 1 and 2 are represented by 
dark, medium and bright areas. 

Looking at each panel in Fig. [3l we can identify cuts for 
which there are broad regions with g = I. For each cut, 
the spin polarization is given by {x\(^\x) — A"^- Figure[4] 
shows these spin components, as a function of e at fixed 
a = .457r (the spins are fully aligned in the z— direction 
for a — .Stt) and as a function of a at fixed e = — I.IJ. 
As one can see, small changes in e (determined by the 
Fermi energy) or in a (determined by the voltage which 
fixes the strength of the Rashba SOI) can cause jumps in 
Sy between large positive and negative values. 



FIG. 4: Spin components Sy and Sz of the polarized propagat- 
ing mode and the evanescent inverse decay length (in units of 
1/L) versus e at a = A5n (left) and versus q/tt at e = — l.lj 
(right). All plots have /3 = .SStt. Data are shown only in 
regions where g = 1. 



In practical situations, one will not use an infinite chain 
of diamonds. As discussed in Refiii, the spin polariza- 
tion is maintained even for a finite chain, provided we 
use adiabatic contacts at the exit. However, this finite 
chain must be long enough so that the evanescent modes 
will decay before the electrons come out. In the sym- 
metric case, cos{qL) is found from solving the quadratic 
equation (A — A)^ — |B|^. For g = 1, one of the so- 
lutions has I cos{qL)\ < 1, and therefore a real q, while 
the other solution has |cos(gZ)| > 1, and therefore an 
imaginary q = in. The bottom panels in Fig. 2] show 
the imaginary part of q for these other modes, denoted 
by K. As one can see, there are broad regions in which 
K > 1, so that a small number of diamonds suffices for the 
evanescent modes to decay. Interestingly, k diverges to 
infinity as a approaches special values, for which the co- 
efficient of coa^ (qL) in the above quadratic equation ap- 
proaches zero. Writing this equation as ax^ +bx + c — 0, 
the solutions for small a are cos(gL) = xi ~ —b/c and 
cosh(KL) = X2 ~ ~b/a. Clearly, k diverges as a ap- 
proaches zero. Using the specific relations which fol- 
low Eq. (O, we identify a = cP + by — b^. It is then 
easy to check that the condition a = is equivalent to 
det W = dot Wt = 0. When this condition is obeyed, W 
and have a vanishing eigenvalue, which means that 
the corresponding eigenvector is completely blocked [see 
Eq. ([3])]. Interestingly, Ref."'^'^ found similar blocking for 
finite chains at special values of the parameters. 
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IV. TIGHT-BINDING VERSUS CONTINUUM 
MODELS 

A. Tight-binding model with many sites 

We now generalize our model, by replacing each edge of 
a diamond by M bonds in series. Consider a single 'large' 
bond, e.g. the bond ab in Fig. [TJ Solving the appropriate 
one-dimensional tight-binding model along this bond, we 
find that the spinor on the site m, ?/'('^^)(m) [0 < m < M, 
with V^^^H*^ = 0) = -iAa and tp^°'^\M) = -0^, is given by 

_ sin[fca(M - m)] . im/A/ , 
^ sin(fcaM) ^^""^ '^'^ 



sin(fcaTO) i(M-m)/j\/ 



(6) 



where C/^^ is the same as in Eq. ([TJ , a = i /Af is the new 
lattice constant, and k is the wave vector for the one- 
dimensional solution, related to the energy e and to the 
new elementary hopping energy J via e = — 2Jcos(fca). 

We next discuss a general node on the network, u. The 
tight-binding equation at this node is 



(7) 



with W = W/7, and with 4A = Ea/j — 4, where now 
7 = ,P /Ef,. The symmetric cases which we described 
above are characterized by the same matrix W. Since 
all the spin physics described above resulted only from 
the matrix W, which does not depend on the energy e, 
all of that discussion will remain unchanged. The only 
effect of adding the internal bonds on each 'large' bond 
appears in the new parameter A, which is now given by 



A = EaEb/{A.f ) - 1. 
As before, the spectrum is determined by 



A = ^± IBI 



(11) 



(12) 



where A — A/^ and B = B/7. For each value of g, Eq. 
p2p determines two values for A, which are independent 
of M and of e. We then solve each of these equations for 
all possible values of fca, and obtain the energies of the 
various bands via e = — 2Jcos(fca). 




where we assumed that all the 'large' bonds are equiva- 
lent to each other, having the same number M of internal 
bonds, the same lattice constant a and the same hopping 
energy J. Substituting Eq. © on the right hand side, it 
is straightforward to rewrite Eq. ([7]) as 



(8) 



within 



j{Nu cos(fcaM) - {Nu - 2) sin(fcaM) cot(A:a) 
e„sin(fcaM)/[Jsin(fca)]), (9) 



where is the number of bonds meeting at site u. The 
generalization to arbitrary bonds is obvious. Equations 
(HI) look exactly like our tight-binding equations ([2]), pro- 
vided we replace e— e„ by E^- Therefore, we might expect 
some similarities in the solutions. 

The equation for in © now has E^ on the left 

hand side, and therefore the dispersionless modes con- 



tain all the solutions of E^, 
e = — 2Jcos(fca), ka = (n 



0. With = 0, this yields 



l/2)n/M. Except for these 
energies, we again eliminate the side site spinors, restrict- 
ing ourselves to the symmetric case, — e^. In this case, 
we also have 7^ = 7^, = 7. It is then convenient to sepa- 
rate the common factor 7 from Eq. and rewrite the 
generalized Eq. ^ in the form 



4AV'„(n) = Wt^Jn ~ 1) + W^Jn + 1), 



(10) 



FIG. 5: Spectra for the diamond chain with AI bonds on each 
'large' bond, with (3 = n/4, (j) = 7r/2, a = An. Top: M = 1 
(left), M = 2 (right). Bottom: M = 4. All energies are in 
units of J. 

For the special case = = = 0, Eq. (fTTj) re- 
duces to A = cos(2fcaM) — .5 sin(2fcoAf ) cot(fca). Ex- 
amples of the spectra for this case and for several val- 
ues of M are shown in Fig. [5l Clearly, the number 
of bands increases with M. However, the basic qual- 
itative shapes of e{q) for all the bands are similar to 
each other. Thus, we expect that within each band 
we would reproduce the filtering properties described 
in Sec. 3. However, the energy scales for each band 
become narrower as M increases. This narrowing is 
most pronounced for the lowest band, which always ap- 
pears below the band of the one-dimensional solutions, 
e < — 2J. Within our tight-binding model, this narrow 
band involves an imaginary value of fc, k = ik, implying 
very small wave functions in the middle of each 'large' 
bond. This imaginary solution results from the solu- 
tion of A = cosh(2KaA/) — .5 sinh(2KaA^) coth(Ka). At 
large M, A becomes negligible, and the two bands con- 
verge to a single band, with coth(Ka) — >• 2 and thus 
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e = — 2Jcosh(Ka) w —2. 3094 J. The bottom panels in 
Fig. [5]show the spectrum for M = 4. The right hand side 
panel zooms on the two lowest bands, demonstrating that 
these bands maintain their qualitative features even as 
M increases. The lowest band does indeed narrow down, 
becoming dispersionless in the limit M oo, where the 
equation A = yields the solution e — — 2cosh(Ka), with 
coth(Ka) = 2. 



B. Continuum model 

In contrast to our calculations, RefsJ^ andi^ used a 
continuum wire model for each 'large' bond on each di- 
amond. Their solution for the spinor at distance x from 
node a on the bond ab is given by 



sin(fca;) 
sin(fcL) 



(13) 



Since the electron on each bond is now free, its energy is 
given by e = h'^k'^ /{2m*), and Bercioux et al. plot e = 
k^. As usual with tight-binding equations, this solution 
coincides with our tight-binding solution ([5]) in the limit 
M oo, keeping a = L/M — » and x = ma. Having 
found these solutions, RefsJ^ andi^ proceed to use the 
Neumann boundary conditions at the nodes: 



E 



a7/;(™)(x) 



dx 



= 0. 



x=0 



(14) 



With these conditions, they end up with equations like 
our ([5]), but with our £"„ replaced by —JNuCOs{kL). 
Thus, they would have A = cos(2fcL), and therefore k — 
(riTTib .5 arccosA)/L. These results, with e — fc^/100, are 
reproduced in the left panel of Fig. [B) Again, the qual- 
itative shape of each band looks similar to ours. How- 
ever, this approach cannot reproduce the lowest band 
with imaginary k which we found for finite M . 




FIG. 6: Spectra for the diamond chain for the continuum 
model (left, with arbitrary energy units) and for the modified 
tight-binding (M = 4, e„ = 2 J, right), with (3 = tx/A, <f) = 
7r/2, a = Att and with energies in units of J. 

As discuss in Ref.— , there is no way to make the tight- 
binding model and the continuum models identical for 
any finite M. However, in the limit M ^ oo one can 
modify the tight-binding model, by setting a site energy 

— 2 J. With this special value, our Eq. (fTT|) becomes 
A — cos(2fcL)-|-sin(2fcL) sin(fca), which becomes identical 
to the expression used in RefSf^i andi^ when a — ^ 0. 
Unfortunately, we know of no good reason to choose this 
particular value for the site energy in the tight-binding 
model. The right hand panel in Fig. [S] shows results of 
the modified tight-binding model for M — 4. Except for 
the dispersionless band below the continuum, we expect 
these results to approach those in the left hand side panel 
as M ^ cx). 



V. DISCUSSION 

ReferenceM already compared the results of a chain of 
diamonds to a single diamond. As also noted in Rcf.^^, 
having more loops in series broadens the parameter re- 
gions which yield full polarization. Reference^^ also dis- 
cussed the conditions for having full filtering on a chain 
of finite length. Here we have extended that discussion 
by showing that there exist broad regions with relatively 
short evanescent decay lengths, so that one can obtain 
filtering with relatively short chains. 

In addition, we have shown that the filtering results 
are quite robust: in addition to the parameters discussed 
in Ref.^'*, the filtering persists upon changing many ad- 
ditional parameters (e.g. the opening angle of each dia- 
mond). We also mention the asymmetry of the spectra 
in the non-symmetric case, which implies that a large 
(positive or negative) bias voltage between the left and 
right hand ends of the device can yield different currents 
in the two directions. 

We hope that the present discussion will stimulate at- 
tempts to realize our filter experimentally. 
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